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Abstract. The goal of the present paper is three- fold. First is to clarify the 
connection between the dominant energy condition and hyperbolicity properties of 
Lagrangian field theories. Second is to provide further analysis on the breakdown of 
hyperbolicity for the Skyrme model, sharpening the results of Crutchfield and Bell and 
comparing against a result of Gibbons, and provide a local well-posedness result for 
the dynamical problem in the Skyrme model. Third is to provide a short summary of 
the framework of regular hyperbolicity of Christodoulou for the relativity community. 
In the process, a general theorem about dominant energy conditions for Lagrangian 
theories of maps is proved, as well as several results concerning hyperbolicity of those 
maps. 



PACS numbers: lllOEf, 1239Dc 
1. Introduction 

In the study of classical field theories, an oft-imposed "physical assumption" is the 
dominant energy condition, which requires that the Einstein-Hilbert stress-energy tensor 
T appearing on the "wood" side [16] of Einstein's equation to have the property that, 
for any future causal vectors X and Y, the contraction T(X, Y) is non-negative. A 
posteriori this condition seems reasonable, in view of the results that can be derived 
from it. Two of the most prominent examples are, of course, the Singularity Theorem of 
Penrose jH] and the Positive Mass Theorem (52j[7T]. However, these results are purely 
results in Lorentzian geometry: that is, Einstein's equation is wholly unnecessary, except 
to transfer the dominant energy condition from the "wood" side to the "marble" side 
of the equation. In other words, those theorems could equally well have been stated 
without reference to general relativity, but merely with certain positivity conditions 
about the Einstein tensor G = Ric — \Rg. 

The typical a priori justification for the imposition of the dominant energy 
condition is some heuristic requiring that the flow of energy be at a speed less than 
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that of gravit)|||, that is, a property on the classical field variously known as finite 
speed of propagation, domain of dependence or causality (see, e.g. [TOJ ) . Perhaps the 
most well-known in this regard is a theorem of Hawking (23] which makes precise this 
notion: that if a matter field satisfies the dominant energy condition and the energy is 
strongly coercive (that is, T(X, Y) — for non-zero vectors X, Y only when the matter 
field vanishes), then if the matter field vanishes in some space-time region, it must also 
vanish in its domain of dependence. Typically, a proof of the domain of dependence 
property uses an energy estimate, which is the curved space-time version of the law 
of energy conservation on Minkowski space-time. And therefore, it is usually expected 
that a domain of dependence property comes hand in hand with the hyperbolicity of 
the associated matter evolution. 

This is, however, not always the case. The principal obstruction is that Hawking's 
theorem only guarantees the finite speed of propagation of "perturbations of vacuum" . 
That is, it essentially only guarantees that the edge of vacuum cannot recede faster 
than the speed of gravity. For semilinear field theories where the kinematics in the 
high frequency limit is always governed by the space-time metric, there can be no 
distinction between perturbations of vacuum and perturbations of a given solution. And 
hence the argument used by Hawking to establish his domain of dependence theorem 
generally indicates hyperbolicity for the associated matter field. For quasilinear field 
theories, however, the strong self-interaction means that the kinematics close to a 
vacuum background can be significantly different from that around a "large" solution. 
And for these types of theories, dominant energy condition is not sufficient to guarantee 
hyperbolicity. 

A prime example of this difficulty is illustrated by two results related to the Skyrme 
model of nuclear physics. Motivated by unstable numerical simulations, Crutchfield and 
Bell demonstrated [H] that highly boosted background solutions to the Skyrme model is 
linearly unstable under perturbations (that there exists an exponentially growing mode). 
On the other hand, it was shown by Gibbons |20j that the Skyrme model in fact enjoys 
the dominant energy condition. In the current paper, a more detailed analysis of the 
breakdown of hyperbolicity for the Skyrme model will be presented. The analysis is 
based on the Christodoulou's regular hyperbolicity framework [TO] . 

(One may also ask whether the reverse implication is true: that hyperbolic matter 
models will always enjoy a dominant energy condition. A simple example showing its 
falsehood is the linear wave equation with a negative potential Du = —u tt + Am = Vu, 
where the potential V < 0. If one were to desire theories without external potentials, one 
can also consider the focusing nonlinear wave equation Du = — \u\ p u, the hyperbolicity, 
or local well-posedness, of which is well-known (see, e.g. [55]). In regions where the \u\ is 
large while approximately constant, the associated Einstein-Hilbert stress-energy tensor 
for either of the above examples violates the dominant energy condition.) 

In the present paper the connections between the dominant energy condition 

| If one is willing to assume a linear or semilinear theory of electromagnetism, then also speed of light. 
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and hyperbolicity will be studied in the context of Lagrangian theory of maps. An 
emphasis will be placed on the insufficiency for mutual implication. The reasons for 
the focus on such matter models are twofold. First is its general applicability. Many 
models of mathematical physics can be cast in the framework of Lagrangian theory 
of maps. Starting from the simple linear wave equation, which can be regarded as a 
map from Minkowski space to the real line or to the complex plane, we can modify 
the target space to a general Riemannian manifold and obtain what is called the 
nonlinear a-model in the physics literature, or the wave-map system in mathematics. 
This system is itself interesting as models in high-energy physics (see [29] and references 
therein) or as symmetry reductions from Einstein's equations in general relativity (see, 
e.g. [HEIHE]). As a semilinear modification to the standard wave equation with a 
geometric interpretation, the well-posedness properties (both global and local, and in 
both small and large data regimes) have been well-studied, see |35 p49U50[l65T - [69"] for some 
recent progress on global regularity and singularity formation in the large data regime, 
and [6l[7l l3^ti3lll5ll[56ll6l] for a sample of classical results in this area. 

In this paper we will focus on quasilinear modifications to the wave-map system, 
whose dynamics is comparatively less well-studied. Such generalisations also have wide 
physical applications, with examples in the nonlinear cr-model hierarchy including the 
Skyrme model [I5j[571|58] , the membrane equation [2UUT], a Born-Infeld type model in 
cosmology [28,38j, and models of hydro- and elastodynamics pll3ll5l[9 l[3Tll5Tl[59ll62] . The 
models of dynamics in a continuous medium are particularly interesting in this context, 
as generally a physical assumption in such models is that the particle world-lines are 
time-like, a condition necessary to guarantee the causality of the matter model. It 
will be shown in this paper that one can construct examples of equations of states for 
which the dominant energy condition is satisfied independently of whether the physical 
constraint is imposed. This further reinforces the idea that the domain of dependence 
theorem of Hawking is only a statement about vacuum perturbations. 

The second reason for considering Lagrangian theories of maps is purely technical. 
For Lagrangian field theories, the Einstein-Hilbert stress-energy can be defined via a 
variational procedure on the Lagrangian density. This allows for general and efficient 
calculations to check the dominant energy condition. Furthermore, in the context of 
Lagrangian theories of maps, the regular hyperbolicity framework of Christodoulou [TO] 
provides a powerful while algebraically simple characterisation of local well-posedness. 
Therefore we will consider only such matter models for ease of discussion. 

The paper is organised as follows. In Section[2]we review the Lagrangian field theory 
of maps and give some examples that have appeared in the literature. In Section [3] a 
geometric method is described for computing the Einstein-Hilbert stress-energy tensor 
for a large class of maps which includes the physically interesting models described 
above. The method provides an easy way to verify the dominant energy condition for 
these maps; we recover the result of Gibbons [20] as a special case. In Section[l]we briefly 
describe the philosophy and method of regular hyperbolicity of Christodoulou [TP] , and 
recall the notion of canonical stress. Here only the basic ideas behind the theory 
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of regular hyperbolicity will be sketched, the focus being on its application. In a 
forthcoming paper with Jared Speck (61], a detailed gentle introduction to regular 
hyperbolicity will be given, along with some simple extensions that were alluded to, 
without proof, in Christodoulou's monograph. And in Section [5] we apply the theory to 
the problem of hyperbolicity of the Skyrme model. 

2. Lagrangian theory of maps 

Throughout we let (M, g) be an m + 1 dimensional Lorentzian manifold, where sign 
convention is taken to be (—,+,+,...); and we let (N,h) be an n dimensional 
Riemannian manifold. M represents the physical space-time (often taken to be 3 space 
and 1 time dimensions, though we make no such restrictions here), while N represents 
the internal structure of the field. In applications, for nonlinear a-models, N is usually 
taken to be a Lie group or a symmetric space; in dynamics of a continuous medium, N 
is the material manifolcijj] with n = m. 

Denote by <ft '■ M — > N a continuously differentiable map. In field theories <ft 
gives the state of the field at a point in space-time, whereas for continuum mechanics <ft 
represents the coordinate transformation between the Eulerian and Lagrangian pictures. 
Then the action of <fi can be used to pull back the metric h onto M as a positive semi- 
definite quadratic form on TM, which we write as 

<f>*h(X, Y) = h(d(f) -X,d<f>- Y) 

where the left hand side is evaluated at a point p G M and the right hand side at the 
point 0(» G N for X, Y G T P M. 

We define the (l,l)-tensor field by composing with the inverse metric g~ x : 

= g- 1 o <p*h . (1) 

We will follow Manton and Sutcliffe [15] and call this the strain tensor for the map, the 
nomenclature taken from the study of dynamics in a continuous medium, where, roughly 
speaking, the trace of on a space-like hypersurface describes the local deformation of 
the material. See p],[3H[62] for more detailed discussions (note that our definition here 
agrees with that of Tahvildar-Zadeh [62] if we impose the physical assumption that the 
map <p admits a space-like simultaneous space, but our definition differs from that of 
Kijowski and Magli [31], who in addition to the above physical assumption, also breaks 
the degeneracy by adding in the square of the particle velocity to make the strain tensor 
positive definite). 

At a fixed point p G M, the tensor defines a linear transformation of the 
tangent space T V M. We can consider its eigenvalues. In the case that g is a Riemannian 

§ In the fluids literature, the material manifold is sometimes only required to be equipped with a volume 
form; but as every paracompact smooth manifold admits a Riemannian metric, and a Riemannian 
metric can realize any volume form through conformal rescaling, it is of no loss in generality to assume 
the material manifold is Riemannian. 
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metric, is a self-adjoint operator on T P M relative to the (positive definite) inner- 
product induced by the metric, and hence all the eigenvalues would be real. For the 
Lorentzian case, the eigenvalues are in general complex. Denote by {Ai, ...,Afc} the 
non-zero eigenvalues, counted with multiplicity. One easily sees that 

k < rank(cZ0) < min(m + 1, n) . (2) 

Recall the elementary symmetric polynomials <7j({Ai, . . . , A^}) given by 

a J ({A 1 ,...,A fc })= J2 ri A ^ ( 3 ) 

l<ai<a2<---<aj<k i=l 

with <j = 1 and <jj = for all j > k. The Oj correspond to the coefficients of 
the characteristic polynomial for D^, with a\ = tr(D^) and a m+ \ = det(D+), and 
are algebraic invariants of the tensor field D*. Furthermore since is a real linear 
transformation, the values of the aj(D^) are all real. By an abuse of notation, we will 
write <jj(D^) when we mean the symmetric polynomials on the eigenvalues of D*. 

In this paper we will only consider Lagrangian field theories for maps <fi : (M , g) — > 
(N, h) where the action integral is given by 

S= [ L(s(<P),a 1 (D^, ( x 2 (D^,...,a rn+1 (D^)dvol g , (4) 

where s : iV — > R is a non-negative scalar function on the internal space N. In particular 
we require the only dependence of the Lagrangian on the field be through the algebraic 
invariants aj(D^) and the value of itself. (Note that by (J2J) any cjj not listed above 
as an argument for L is not dynamical.) Physically the dependence on aj(D^) and 
not other components of corresponds to the assumption that the laws of physics 
are locally invariant under a Lorentzian rotation in 0(1, m) of the tangent space T V M 
that fixes the kernel of D^. In relativistic elasticity where a space-like simultaneous 
space is imposed, the relevant subgroup of 0(1, m) to consider is the orthogonal group 
0(n), and this condition on the action corresponds to the assumption that the material 
is homogeneous, isotropic, and perfectly elastic [62]. The function s corresponds to the 
entropy per particle in relativistic elasticity, and plays the role of the symmetry-breaking 
mass term in the Skyrme model. 

Observe that since the action itself only depends on D®, which is defined through 
only the first derivative of the map 0, the equations of motion given by applying the 
Euler-Lagrange equations to (0J will be quasilinear second-order partial differential 
equations. 

Some explicit examples of Lagrangian field theories from the literature that fall in 
this class include: 



Wave maps can be described by the Lagrangian function 
L = a 1 ( J D*) = |V0| 2 
where the norm is taken relative to both g and h. 



(5) 
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Relativistic elasticity typically makes the additional assumption that d<p is onto, 
and its kernel time-like (the space-like simultaneous space assumption). Then 
interpreting s as the entropy per particle, the action (jlj) is the general form for 
studying a homogeneous, isotropic, and perfectly elastic solid. The case where 
L = L(s,a n ) is the special case representing relativistic fluids, a n = \4>*dvolh\ 2 
being equal to the squared g-norm of the pullback of the volume form on N. 

Skyrme model classically assumes (M,g) to be the 4-dimensional Minkowski space 
with (N, h) the standard 3-sphere, though we can also consider the model without 
such constraints. The Lagrangian function is (up to rescaling) [571 [58] ( see a l so 

MM) 

L = <t 1 (D*) + <t 2 (B*)+s(<P) , (6) 

where <72(-D^) = \d<f> A d(p\ 2 , and s(0) = mdisth((j), 4>o) is a multiple of the geodesic 
distance of from a fixed point 0o £ N, the constant m being the mass parameter. 
In general, Lagrangians of the form L — a\ + . . . are the nonlinear cx-models. 

Born-Infeld type models fix a large constant b > 0. Restricting the allowed maps 
to those whose eigenvalues of have real parts at least —b, the action is given by 

L = y/det(b ■ Id + W) - v/det(6 ■ Id) . (7) 

This is the cx-model analogue of the Maxwell-Born-Infeld model of nonlinear 
electrodynamics, see also [20] . 

Membranes as considered in [SHE] can be described by setting (M, g) to Minkowski 
space of some fixed dimension, and (N, h) to the real line with canonical metric, 
and taking 

L = y/l + a^D*) (8) 

while restricting consideration to those <p f° r which L is well-defined. It can also 
be viewed as an analogue to the minimal surface equation for embedding into 
Minkowski space. 

3. Dominant energy condition 

In this section we describe some sufficient (but far from necessary) conditions on the 
Lagrangian function L of the form considered in the previous section that guarantees 
the associated Einstein-Hilbert stress-energy tensor satisfies the dominant energy 
condition. In isolated cases (fluids and elasticity with the space-like simultaneous surface 
assumption, or exactly the classical Skyrme and Born-Infeld models (20]) the results 
are known before. The novel contribution in this paper is given below in Proposition 
13 It can be used in a unified algebraic framework applicable to all theory of maps 
described by an action of the form (TJ]), making unnecessary the ad hoc computations 
through explicitly evaluating the eigenvalues used in e.g. [20]. It is worth remarking 
that those types of computations rely on a genericity argument to diagonalize a positive 
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semidefinite quadratic form relative to a Minkowski metric, a procedure which cannot be 
carried out when admits a null eigenvector. By working geometrically and tensorially, 
the computations described below avoids this difficulty 

We start by reviewing some definitions. Recall that the Einstein-Hilbert stress- 
energy tensor T G r(T 2 °M) for a Lagrangian field theory is given by a variational 
derivative for the Lagrangian density relative to the inverse metric, 

T^d^:=2Mpl=(^-iL 9 )^diM. (9) 

Definition 1. The stress-energy tensor T is said to obey the dominant energy condition 
at a point p G M if WX G T p M such that g(X, X) < 0, the following two conditions are 
satisfied 

T(X,X)>0 (10a) 
[Tog- 1 oT](X,X) < (106) 
unless T vanishes identically. 

Remark 2. The definition is equivalent to the classical statements (see, e.g. section 4.3 
in [23] or chapter 9 of [70] ) of the dominant energy condition. Observe that (11061) gives 
that the vector g^ 1 oT oX is a causal vector for any time-like vector X, and fllQqj) gives 
that the vector g^ 1 o T o X has opposite time-orientation as the time-like vector X. 

The set of future-pointing time-like vectors form a convex cone; hence we have the 
following technical lemma, applicable to all Lagrangian field theories, not just those 
described in Section [2j 

Lemma 3. Let F = F(x\, . . . , Xk) be a continuously differentiable function of k real 
variables. Assume that F is concave, F(0) > ; and that diF > for each 1 < % < k. 
Let Li, 1 < i < k denote a collection of Lagrangian functions, and let Ti denote their 
corresponding stress- energy tensors. Suppose Tj each separately obeys the dominant 
energy condition, or, equivalently, the vectors Y\ = g^ 1 o Ti o X are all past-causal for 
any fixed future time-like X . Then the stress- energy tensor T for the Lagrangian formed 
by L = F(L\, . . . , Lfe) also obeys the dominant energy condition. 

Proof. The stress-energy tensor T can be written, using (Q, as 

i=l y i=l i=l 

Now considering g~ 1 oToX , the first term in the above expression contributes 9iF-Yi. 
By assumption, this is a convex combination of past-causal vectors, and hence is past- 
causal. For the second term, since g^ 1 o g o X = X, to show that it is also past-causal 
it suffices to show that 

k 

F>J2 9 i F ' L i ■ 

»=i 

But this follows from the fact that F is concave and F(0) > 0. □ 
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Remark 4. That F is required to have nonnegative partial derivatives represents the fact 
that each of the Lj's contribute nonnegatively to the energy. That F(0) > states that 
there is no negative vacuum energy. Both conditions are therefore natural and necessary 
for the total Lagrangian L to have positive energy density, if the Lj's are taken to be 
independent. The concavity condition is technical. It appears naturally in the proof, 
but can potentially be relaxed if more is assumed on the individual Lj's. 

We will apply Lemma |3] to the following proposition, which is the main 
computational result of this section. Observe also that for L = s(<j>), its corresponding 
stress-energy tensor is T — — |s • g, and by the assumption on the positivity of s, obeys 
the dominant energy condition. 

Proposition 5. For L = o-j(D^), T obeys the dominant energy condition. Furthermore, 
T = at a point p if and only if j > rank(<i0| p ). 

Before giving the proof, we need to review some linear algebra. Consider a real 
vector space V. Let A be a linear transformation on V. Then A naturally extends to a 
linear transformation, which we denote A^ , on A 3 (V), the space of alternating j-vectors 
over V. A classical result in linear algebra is that o~j(A) is proportional to tr^irv) ^ ■ 
Now, letting V = T p M and A = = g^ 1 o <fi*h, we observe that 

or, to put it in words, (D^)^ is obtained from first taking the induced metric hP on 
alternating j- vectors in T^^N, pulling it back via 0, and composing it with the induced 
metric (g~ 1 )^ for the alternating j-forms. In index notation, this can be written as 

[{D^t;:% = g blC1 ■ ■■g b ^{^h) ai[cl \{<p*h) a2 \ C2 \ ■ ■ ■ m) aj _ x \ Cj _^h) aj \ Ci] 

where the bracket notation in the indices denotes full anti-symmetrization of the 
{ci, . . . , Cj} indices. For a Lagrangian proportional to aj, we can assume 

L = [{D*ff a x;% = s ai[ci1 • ■ ■ g a ^M*h) aiCl ■ ■ ■ {4>*h) ajCj . (li) 

It is simple to check, using (D^) = diag(—l, 1, 1, . . .) that the above expression has the 
correct sign: that L defined thus is a positive multiple of aj. 

One can also arrive at fflTT) purely from a linear algebra point of view. Let pj be 
the power sum 

k 

P ,({\ 1 ,...,x k }) = Y J K ■ 

i=l 

Recall Newton's identity 

j 

3 ■ <Tj = ^(-ly^CTj-iPi 
i=l 

which allows us to express aj as a rational polynomial in pj's. Now, by definition, it is 
clear that 

Pj {D*)=tr[{D*Y] 
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i = <r i *'''!p*K:ZW 



h) aiCl ■ ■ ■ (4>*h)c 



Newton's identity reduces to a generating condition for E based on the Kronecker 5 
symbols, 



E c b = Si , 



Cl ... Cj 
J Cj b 1 ...b j 



E, i \i— 1 rp c i--- c j-i K c j-i r c j-*+i r c i-i+i 



i=l 



a...cj 
'j 



A direct computation which we omit here shows that then in fact the invariant E bi "' b 
is a positive rational multiple of the generalized Kronecker symbol 8^ '"^ , from which 
we recover ffTTl). 



Proof of Proposition^ We need to show that g^ 1 o T o X is past causal for any future 
time-like vector X. Since T is tensorial, we can assume X has unit length. Using the 
expression for (a positive scalar multiple of) Oj given in (iTTj) . we can write T(X, •) for 
L = <jj in index notation: 

T ab X b = IV N ■ • • 9 aAcj] {<t>*h) ab ■ ■ ■ {cj>*h) a . Cj - l -g ab X h L (12) 

Take an orthonormal basis for T p M relative to g. Since we assumed X unit, let 
eo = X and {ej}i<j< m all space-like. We can take j < m + 1 as otherwise T is identically 
0. Then we notice that a basis for hP(T p M) is given by 

{e A e ai A • ■ ■ A e aj _J 

We write the first set as and the second set as A|. Using the normalization that 
v /\w = v ® w — w ®v, we find that each of the element in A J ± has norm — j\ while the 
elements in A| has norm j ! . 

To show that T(X,X) > 0, we observe that under the expansion ffl2l . the first 
term corresponds to 

while the second term corresponds to 

/ \ 

So summing them gives 

/ \ 



v 



/ 
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which is non-negative by the fact that <p*(h^) is a positive semi-definite quadratic form 
on A J (T p M). Furthermore, observe that since A| U Aj_ is a basis, its push- forward 
0*Ajj U 0*A^_ spans A J '{(p^TpM) C A J (T^N). Thus by the fact that ft, (and hence the 
induced metric is positive definite, we conclude that when T(X,X) = 0, necessarily 
A J (0*T p M) = {0}. This proves the assertion that T vanishes only when j > rank(c?0). 
To show (llO&p . we observe that 

m 

X a T ac rt b X b = -T(X,X) 2 + ^T(X,e,) 2 . 

i=l 

The first thing to note is that T(X, e,) does not have any contribution from the second 
term in (fl2i) . For the first term, a quick computation shows that T(X,ei) corresponds 
to 

<f>*{h ij ){e Av,eiAri) 

SO 

m 

|^r(x, ei ) 2 |<(^iT(x, ei )|) 2 
i=i 

m 

<(E E I^^KeoA^e.A^I) 2 

i=1 »7gA|' | ~ 1 ,eiAr?^0 
1 m 

1 m 

= i( E E^')( e * A ^ A ^)) 2 

= T(X,X) 2 

and therefore (1106|) is satisfied. □ 

As an immediate application of Lemma [3] and Proposition [5j we have that the 
Skyrme model and the Born-Infeld model described in Section [2] obey the dominant 
energy condition: it suffices to check that -F Skyrme (s, 0"i, <t 2 ) = s + cr 1 + cr 2 and 

Fbi(o"i, . . . , <J m +i) = \Jj2T +1 b m+1 ~i(jj — Vb m+1 are concave, satisfies F(0) > 0, and 
has positive partial derivatives, conditions which are easily seen to hold. Therefore we 
recover the result of [20] in dimension m = n = 3, and also extend it to arbitrary 
dimensions m,n. In fact, we have 

Theorem 6. Given a Lagrangian theory of maps with action given by a sufficient 
condition for its Einstein- Hilbert stress-energy tensor to obey the dominant energy 
condition is that the Lagrangian function L in be continuously differentiate with 
nonnegative partial derivatives on its arguments, be concave, and satisfy L(0) > 0. 
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Remark 7. We stress again here that our notion of Lagrangian theory of maps, as 
discussed in Section (2J is very general, and includes the general forms of Lagrangians 
used in relativistic elasticity and relativistic fluids. Therefore it is remarkable that for 
these theory of maps, the assumption of a time-like particle world-line is inconsequential 
insofar as the dominant energy condition is concerned: in particular a tachyonic fluid 
will still obey the dominant energy condition, in direct contradiction to the intuition 
often presented as the justification for the dominant energy condition. For a concrete 
example, consider the fluid Lagrangian with some large constant b on Minkowski space, 

L = y/b + a 3 (D<l>) . 

By Theorem [6] it obeys the dominant energy condition; locally, however, the "fluid" 
moving with a constant "velocity" larger than that of gravity is a solution, e.g. 
(f) : ]R 1+3 — y R 3 where (fi(t, x, y, z) = (t, x, y) is a solution to the equations of motion. As 
we will see in Section |5] for the Skyrme model, however, there is an instability (i.e. a 
violation of hyperbolicity) associated to these tachyonic matter. 

4. Regular hyperbolicity and the canonical stress 

The question of hyperbolicity for a system of partial differential equations is generally 
synonymous with whether the system admits a locally well-posed Cauchy problem for all 
smooth initial data. In the context of linear systems of constant coefficients, it is a well- 
known theorem [251 Ch. 5] that a necessary and sufficient condition for hyperbolicity 
of a system P(d)<f) = is the hyperbolicity for its polynomial symbol P(£). (That the 
Cauchy problem is well-posed for all smooth initial data is crucial for the necessity. 
See Remark [8] below.) Two difficulties arise when trying to apply this theorem to 
non-linear or variable coefficient systems. First is the general difficulty of checking 
the hyperbolicity for a symbol, which requires computing the singular locus of P(£) _1 . 
Second is the technical difficulty that concrete quantitative estimates are required for 
iteration arguments used in these type of problems. 

Many stronger notions of hyperbolicity exist in the mathematical literature (see 
|12[|3T] for some samplers); the common theme to all is that they provide sufficient 
conditions for local well-posedness of the non-characteristic Cauchy problem. In general, 
however, these conditions are not necessary: that Leray hyperbolicity |39j is a subset 
of symmetric hyperbolicity [18] is well-known; that the Maxwell-Born-Infeld model 
can only be seen as symmetric hyperbolic after an augmentation procedure [2| [53] 
strongly suggests that the notion of symmetric hyperbolicity does not directly capture 
all hyperbolic systems. (It is perhaps interesting to note that the hyperbolicity of 
Maxwell-Born-Infeld system can be directly treated within the regular hyperbolicity 
framework |60j.) In this section we will describe the regular hyperbolicity framework, 
which provides another sufficient condition for local well-posednes^jj]. To guarantee local 

|| That this framework is not better known in the community is perhaps due to the dense mathematical 
language in Christodoulou's monograph |10) . The author hopes to provide here, by way of a summary 
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well-posedness, the general technique common to these methods is that of the L 2 energy 
estimate. 

4-1. Energy estimates 

The energy estimates that we will need are estimates of the following form: assume 
for the current exposition that M is diffeomorphic to M 1+m , and let H k denote the 
L 2 Sobolev space with k derivatives on constant time slices diffeomorphic to R m . The 
energy estimates are inequalities 

WWh-V) < ||0||^(O)e c * VO < t < T,C = C(T, \\<P\\ H k(0)) . (13) 

(That we are only interested in such L 2 energy estimates is because we look for estimates 
generally applicable to hyperbolic equations. In particular, we require the estimates to 
hold for the linear wave equation. By a theorem of Littman [42] , inequalities of the form 
Q13p for the wave equation can only hold in L 2 based function spaces.) Inequalities of 
the form ( {TBI guarantee local existence up to time T. 

To obtain such energy estimates, the usual method is via the divergence theorem. 
Let J = J(t, x, 4>, ... , d k <p) be a vector field, depending on up to k derivatives of (f> which 
we assume solves some partial differential equation. The divergence theorem states that, 
for n = d t the normal vector field to the slices 

— / J(0, x) ■ ndx + / J(t,x)-ndx= / / div(J)dxdt . 



The inequality ({TBI) would hold after one application of Gronwall's lemma, if we can 
guarantee the following 



C \\<t>\\ H * (t) < J J(t,x) -ndx<C\\(j>\\ H k(t) forsomeC>l, (14a) 

J div{J){t)dx < C'\\<j>\\ H k{t) for some C > . (146) 

A J verifying ( 11 4 aft and ( 1146ft will be called a compatible energy current. (In [10J, 
Christodoulou use the term "compatible current" to refer to any J verifying (114fo|) . 



The use of the word "energy" here is meant to reflect the imposition of the additional 
positivity condition (II 4 aft .) 

For Ijl4l to hold, it is necessary that we apply the equation: by the chain rule, 



the divergence of the vector field J which depends on k derivatives of will depend on 
k + 1 derivatives of 0; it will be absurd to be able to control its integral by something 
depending on fewer derivatives. But by suitably applying the equations of motion for 
0, we can convert top-order derivatives to lower-order ones and satisfy the inequality. 

of the main ideas behind the theory, an advertisement for this simple yet powerful technique that is 
well- adapted for use in Lagrangian field theories. In a forthcoming paper [61], Jared Speck and the 
author will give a more accurate and detailed introduction to regular hyperbolicity and fill in some 
material omitted in 101. 
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4-2. Regular hyperbolicity 

Christodoulou provided in his monograph [10J a robust and geometric way for obtaining 
compatible energy currents for any given Lagrangian theory of maps, using techniques 
similar to those used by Hughes, Kato and Marsden in a non-geometric context |26j. 
In this subsection, we quickly review the key points of the theory. The main results 
obtained by Christodoulou are that 

(i) Vector fields satisfying (I14frj) . up to a total divergence term, and up to Noether 
currents coming from symmetries of the target manifold, generically arise from 
contractions of arbitrary vector fields against what is called the canonical stress 
tensor, which can be obtained algorithmically from the equation of motion. 

(ii) The existence of vector fields satisfying (I14q[) depends only on the causal properties 
of the canonical stress tensor. 

Let us begin by summarising the construction of the canonical stress tensor. Let ip 
be a function taking values in W 1 . We will use index notation where components of R" 
are indicated by capital Latin letters. Assume ip solves a system of second order partial 
differential equations of the form 

m%d 2 ab ip B = F A (iP,diP) (15) 

where lowercase Latin letters denote indices on the space-time manifold, and m a A b B is 
some field of coefficients, symmetric in a, b, and symmetric in A, B. Consider the tensor 
field defined by 

Z%\ c d := m%5 c d - m%5 a d - m A c B 5 b d , (16) 
it is a direct computation to show that the tensor field 

QMS := -Zf B \ d d a i> A d b r (17) 
has the property 

^QMd = -(5 c ^|S)^^ + 2(5 d B )F B ^,^) • (18) 

Therefore, for any vector field X, the vector field Q c d X d has the property that it depends 
on 1 derivative of ip, and so does its divergence. This tensor is the canonical stress 

tensor associated to the solution ip of the equation (fl5l) . Notice it only depends on dip 
and the coefficients for the highest order derivative term, m AB 

In general, for a Lagrangian theory of maps described in Section [2j the Euler- 
Lagrange equation will take the quasilinear form 

m%(x, 0, ,90) V> B = F A (x, 0, 90) . 

By commuting in further derivatives, we see that the partial derivatives ip A = V a A 
all satisfy equations of the form (fT5l) with the same m a A b B . Therefore, with judicious 
applications of Sobolev embedding theorem, we see that for any fixed X, the vector 
field 

J c :=X c |0| 2 + ^Q[V a 0]^ rf 

|a|<fe 
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satisfies (1146|) . Note that thus far the construction of J is an algebraic and algorithmic 



statement about Lagrangian theory of maps. 

To obtain hyperbolicity, it is necessary to also satisfy ( 114 aft . It was shown by 
Christodoulou that given a foliation S 4 of M by level sets of a function t, that a sufficient 
condition for J (defined via the vector field X) to satisfy (114 aft is for the following to 
hold 

(i) m^ B (dt) a (dt)b is a negative definite matrix; 

(ii) m^g^a^b is a positive definite matrix for all non-zero one-forms satisfying = 0. 

(Please refer to pU] for proof.) We shall call functions t that satisfy the first property 
time functions and vector fields X that satisfy the second property observer fields. 
For semilinear equations where m°^ B = g ab tiAB, these notions agree with the that of 
time functions and timelike observer fields in general relativity. In the quasilinear 
case, the time functions and observer fields form a replacement for the usual causal 
structure of the Lorentzian metric for governing the kinematics in perturbative analysis 
of solutions. One can also analogously define the notion of global hyperbolicity, domain 
of dependence, and maximal development of initial data relative to this replacement 
causal structure. (This formulation is, in particular, used in Christodoulou's seminal 
work on shock instabilities of the Euler equation [11].) In [10] it is claimed without proof 
that the existence of a time function and an observer field is sufficient to guarantee the 
local well-posedness of the Cauchy problem with data prescribed on a level surface of 
t, a proof will be supplied in [6]]; for the non-geometric scenario working over a fixed 
coordinate system in Minkowski space, with X — d t , a proof is available in [26J. 



4-3. Breakdowns of hyperbolicity 

When the existence of time functions and observer fields fail, the regular hyperbolicity of 
the system breaks down. This in particular implies the non-existence of general energy 
estimates by the work of Christodoulou [10], and hence the impossibility of applying 
the usual iteration method to obtain local existence and uniqueness of solutions to the 
Cauchy problem. While this does not prove the lack of hyperbolicity for those systems, 
here we provide some heuristic arguments as to why the lack of regular hyperbolicity is 
indicative of a lack of local well-posednes^J]]. 

To illustrate the different modes of breakdown, we first consider the linear equations 
given by 

eo«tt + eiUn + £2^22 + e 3 «33 

for some scalar a on R 4 , with e* taking values ±1. The coefficients m°^ B for this equation 
take values in 1 x 1 matrices, i.e. scalars, and we have that 

m ab = diag(e ,e 1 ,e 2 ,e 3 ) . 

% We should recall at this stage that we consider well-posedness in the sense of Hadamard [21] : the 
existence of suitably regular solutions, the uniqueness of said solutions, and the continuous dependence 
of the solutions on given initial data. 
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The usual case of the wave equation is given by — e = t\ = e 2 = e 3 = 1. It is clear that 
then for a given covector 

f = £ dt + ixdx 1 + i 2 dx 2 + £ 3 dx 3 , 

the requirement for m ab ^ a ^b to be negative definite (that is, negative as a scalar) is the 
usual condition that £ is time-like with respect to the metric m ab : £q > £f + £| + £|. 
Similarly, for m ab £ a £& to be negative, we just reverse the preceding inequality. Hence the 
wave equation is hyperbolic if we choose t to be a time-like (in the usual sense) foliation 
and X to be any time-like (again in the usual sense) vector. 

The standard example in which one cannot construct any time-function is the case 
of Laplace's equation, where e# = 1, and for any covector m ab ^ a ^b is positive. For 
elliptic problems, it is well-known j22j[63] that the Cauchy problem is ill-posed: there 
cannot be continuous dependence on initial data. By comparison, we shall say that 
a Lagrangian theory of maps have an elliptic type breakdown of hyperbolicity if one 
cannot construct any time-function even locally. Note that elliptic type breakdowns 
are not the same as the equations forming a bona fide regularly elliptic system, which 
requires the "Legendre-Hadamard condition" (TO] that for all covectors £, m a ^ B ^ a ^b be 
positive definite matrices. What we call elliptic type break down only suffices that none 
of those matrices be negative definite, and in particular mixed signatures will imply 
break down. 

The case where one can construct a time-function but not any observer fields we 
refer to as ultrahyperbolic type breakdown; the name is taken from the canonical example 
of the ultrahyperbolic equation 

-U U ~ «11 + «22 + U33 = 

where any foliation with normal covector £ satisfying £q + £f > £| + £f contributes a 
time-function, but the trace of m ab to any three-plane is indefinite. The ultrahyperbolic 
equation have infinite speeds of propagation [27J, and can be checked by the theorem 
alluded to in the beginning of this section to have a non-hyperbolic polynomial symbol, 
and hence cannot admit well-posed Cauchy problems. (The instability of ultrahyperbolic 
equations has also been considered on physical grounds in the literature, see e.g. [T5].) 
Remark 8. The ultrahyperbolic equation illustrates an important connection between 
hyperbolicity and finite speed of propagation. As mentioned, the ultrahyperbolic 
equation does not admit well-posed Cauchy problem for all smooth initial data. It was 
however pointed out by Craig and Weinstein [13] that, if one were willing to impose "non- 
local constraints" (in their case a correlation on the admissible space-time frequencies 
of the waves), the Cauchy problem can be well-posed in the restricted class. One should 
think of the non-local constraints as circumventing the instabilities caused by infinite 
speeds of propagation. For semilinear equations, it may be possible to impose such 
constraints a priori and globally (provided these additional constraints are compatible 
with the equations of motions; the constraint given in [13, §2] is not preserved if one were 
to modify the ultrahyperbolic equation by a power nonlinearity |m| p m); the situation for 
quasilinear equations is much less clear. 
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4-4- Applying to Lagrangian theory of maps 

A direct computation shows that the coefficient tensor m°^ B can be obtained as the 
second variational derivative of the Lagrangian function relative to the field velocity. 
That is: 

m%(x, 0, d<f>) = 6da(p f sdb(pB L^ 0, dd>) . (19) 

Noting that the positive (negative) definite matrices form a convex cone inside the space 
of matrices, we see that if t and X are simultaneously time functions and observer fields 
for a collection of Lagrangian functions L iy then they will also form a pair of time 
function and observer field for any convex sum of the Lj's. 

Turning now to the individual cases L = aj(D^) as described in Section 121 we can 
again use (TTTj) to facilitate computations. Directly we get 

• (■t>'h)« 1 _„ s _ l Slfi l h BD + 2(j - l)4 J _ 1 9„ J _ l) ^S(a^' 4 '^,^sftA'D 

Now consider an arbitrary time-like unit covector £, and let g be the restriction of g to 
the orthogonal complement of £ (in particular it is a positive semidefinite bilinear form), 
we see that we can write rn b BD ^d in the following way 

mgotei = -2g a ^ ■ ■■g a >- llc i- l] (<P*h) aiCl ■ ■ ■ (<P*h) a ._ 2Cj _ 2 

■ [(fh^^hBD - d c ^4> A d aj _ l( p c h AB h CD ] (20) 

after carefully counting the various permutations of indices leaving the expression non- 
zero. Using a Cauchy inequality on the term now in the brackets, we see that it is positive 
semi-definite; it will be positive definite whenever d(f)\ pe M '■ T P M D £- L — > T^N is has 
rank at least 2 (that is, the tangent space map d(p restricted to the vectors in the kernel 
of £ has rank at least 2). Therefore, we have the following theorem. 

Theorem 9. Let the Lagrangian L be a convex linear combination of aj(D^) and s(4>). 
And let t be an arbitrary time function relative to the underlying metric g a \>. Then 

(i) if the coefficient for o~\ is non-zero, then t is a time function for the equations of 
motion; 

(ii) if d(f) restricted to the level sets oft has rank > max(j — 1,2), then t is a time 
function for the equations of motion; 

(Hi) failing both of the above, t is borderline degenerate, but the corresponding 
m AB(dt) a {dt)b is positive semi- definite. 

In particular, L cannot have a bona fide elliptic type breakdown of hyperbolicity. 

Noting that the domain of dependence in theory of hyperbolic equations is 
determined precisely by the admissible time-functions [ffilEH], we have the following 
interesting corollary, which implies that a Lagrangian theory of maps with the dominant 
energy condition (compare Theorem [6]), when in fact hyperbolic, cannot propagate 
perturbations faster than the speed of gravity. 
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Corallary 10. Let L = L(s, <7i, . . . , aj) be a Lagrangian function with the property that 
L is non-decreasing and concave in its arguments. Assume further that for some k, 
the partial derivative of L with respect to at is bounded from below, and that either (a) 
k = 1 (b) admissible solutions have rank(d<p) > max(A;,3). Supposing the dynamics 
derived from L is regularly hyperbolic, then the domain of dependence for solutions must 
be strictly contained within that of the linear wave equation on the space-time (M, g) . 

Proof. Using the energy estimate, it suffices to show that any time function t of the 
underlying space-time is a time function for the Lagrangian theory (see pJJ §5.3] for a 
discussion). Using ( 1T9"1) we write 

nh d 2 L 5ai 5oj dL 5 2 aj 
m% - 



doidoj 5d a (j) A 5db4> B daj 5d a <t> A 5db4> B 
After contracting with (dt) a (dt)b, the first term is positive semi-definite using the 
concavity of the Lagrangian functional. The second term is a convex sum of positive 
semi-definite matrices, and by Theorem |9] and the hypotheses, at least one of the terms is 
positive definite. Hence t is a time function for the Lagrangian theory of maps described 
by L. □ 

Remark 11. We can compare the above corollary to the situation in isentropic fluids. 
After choosing an equation of state, one can write the pressure p of the fluid as a function 
of the proper energy density p. In this regime the well-known criterion for hyperbolicity 
is that the speed of sound \/ dpjdp is positive and real, while the criterion for the fluid 
flow to be causal is that the speed of sound is less than the speed of gravity (which we 
can choose to equal 1) [62] ■ The proper energy density is in fact the Lagrangian function 
for a fluid. Writing r = (a n (D^))- 1 / 2 for the volume per particle, the pressure can be 
given [10] by p = —4-{rp). Then with a bit of elementary calculus, one sees that the 
causality conditions are equivalent to the following statement on L = L(a n (D^)): 

> 2a n L" > —L' . 

The first inequality guarantees the speed of sound is smaller than that of gravity, and 
it states precisely that L is concave, agreeing with the corollary above. The second 
inequality guarantees hyperbolicity; compare with the discussion below. Note that a 
necessary condition for the second inequality to hold is that L is increasing as a function 
of (T3, which is also a part of the hypothesis in the above corollary. 

The analogous statement for observer fields, however, does not hold. Letting r\ now 
be an arbitrary space-like unit covector, and let X be a unit time-like vector orthogonal 
to rj. Denote now by g the restriction of g to the orthogonal complement of rj and the 
metric dual of X. We see that the analogous statement to ( |20l is the following 

m b i Md = z~g ai[ci1 ■ ■■g aj - llCj - l] {<t>*h) aicl ■ ■ ■ (<P*h) aj _ 2Cj _ 2 



[(<P*ti) aj _ l0j _ih B D ~ Qcj-i^daj-ilp hABhcD] 



-^■■■r jMci -H<t>*h) aicl ---{<t>*h) a 



j- 2 Cj-2 
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• [x<p A x<p c h AC h BD - X$ A X<fh AB h CD \ 

zi r i cl i... rJ - 2 | Cj - 2 ] (0 */ l)aici ... ( ^/, )aj 

X<p E X<f) F h EF [(<t>*h) a ._ 2C ._ 2 h BD - d c ._J J 



' A d aj _ 2 (jfh AB hcD\ 



(21) 



where the first term on the right hand side is again positive semi-definite by Cauchy's 
inequality, while the second and third terms are a priori negative semi-definite. In 
particular, comparing the first and third terms, if is large compared to g ac {4>*h) ac , 

the bilinear form m b BD ribi r ]d can easily have negative eigenvalues; one can compare this 
to the conclusion drawn in p3] about instabilities of the Skyrme model. We are forced 
to conclude that a general Lagrangian theory of maps is susceptible to ultrahyperbolic 
type breakdowns of hyperbolicity. 

On the other hand, notice that in the case Xcf) = 0, the form m BD r}i)Tid is indeed 
positive semi-definite, and furthermore positive definite if d<p satisfies an analogous rank 
condition as in Theorem [9j Also note that in the case j ; = 1 (the semi-linear case), 
positive definiteness of m b BD r]brid holds without need of any assumptions on X<p, as the 
problematic terms do not appear. 

Remark 12. As seen in the discussions above, the term <j\ in the Lagrangian always 
introduce a factor that is regularly hyperbolic. This in fact has a stabilizing effect on 
the hyperbolicity of the field theory. An example will be given in Section El where 
hyperbolicity of the Skyrme model persists into a regime where the particle velocity 
(using a fluids interpretation) exceeds that of speed of gravity. 

4-5. Canonical stress versus Einstein- Hilbert stress-energy 

Thus far we have seen that for o~j, the dominant energy condition must hold; on the 
other hand, the computations above show that the analogous statement for the canonical 
stress need not hold. That is: for perturbations ip over a solution cf) to the equations of 
motions determined by aj(D^), the canonical stress need not have the property 

"QVltfdX 6c > for any time-like vector X and time-like covector £ with < 0". 

However, it is also clear from the definitions that the canonical stress tensor and the 
Einstein-Hilbert stress-energy tensor in fact agree in the highest order derivative terms 
for the semilinear case L = ai(D^), in which the equations of motion are always regularly 
hyperbolic. What is, then, the difference between the two stress tensors? 
In fact we have the following: 

Proposition 13. For L = aj(D^), the canonical stress tensor Q[<p] corresponding to 
the solution itself agrees with g~ x o T , the Einstein-Hilbert stress-energy tensor. 

The proof follows from a direction computation of Q[<f)] d X d using the formulae given 
in ( fl9l) et seq., and comparison against ( fl2l) . We omit the details here. This proposition 
shows that for Lagrangians composed of convex linear combinations of the o~/s, the 
dominant energy condition is sufficient to guarantee the existence of a compatible energy 
current that controls the solution itself, but not its higher derivatives. This is, of course, 
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unsurprising, as the Einstein-Hilbert stress-energy tensor is divergence free by definition, 
and hence for any time-like vector field X, g^ 1 o T o X is a compatible current. 

It is tempting to propose that, in view of Proposition [131 regular hyperbolicity 
implies the dominant energy condition. This is, however, not necessarily the case 
in general. Returning to (1 14 aft of the fundamental energy estimate, we see that it 
only suffices that the compatible current is comparable with the Sobolev norms in 
an integrated sense, and indeed, the conditions on regular hyperbolicity encoded in 
positivity and negativity conditions on m AB only guarantees that much. By way of an 
example, we consider a linear equation on M 1+2 with a two-dimensional target manifold. 
Consider the following set of matrices 

(22a) 
(225) 
(22 c) 



(22d) 



(22 e) 

and let m AB = m AB — e5 ab SAB for some e« 1. Now, by construction, m A ° B is clearly 
negative definite, so the usual function t is a time function. A direct computation shows 
that rh A b B r] a r)b is positive definite for any non-zero rj that satisfies r](dt) = 0. Hence for 
sufficiently small e, the coefficient matrix m AB is regularly hyperbolic, and the equation 

<B^ B = 

has a well-posed initial value problem for finite energy initial data. However, the 
canonical stress tensor (J^lo is not point-wise positive definite for all initial data! Take 
an arbitrary smooth initial data with compact spatial support on t = such that inside 
the ball of radius 1, the data takes the following explicit values: 

ij\0,x,y) = y , (23a) 
?P 2 (0,x,y) = -x , (236) 
d^ A (0,x,y) = . (23c) 

It is easy to check that for this initial data, ih a AB d a ip A dbip B = in a spatial neighborhood 
of 0. And hence in that neighborhood the energy density for Q[^]o is i n f ac t point- wise 
negative. 

This lack of point- wise positivity is a reflection of the internal structure of the target 
manifold N . Indeed, if N were one-dimensional or effectively one-dimensional (i.e. the 
coefficient matrix m AB = g ab h A B is separable), two waves that are spatially coincident 
and travels in the same direction must interact, and hence a point- wise measurement 
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can capture correct notion of energy. But when N has an internal structure, two 
waves that are spatially coincident and travels in the same direction does not have 
to interact strongly, since they maybe orthogonal on N. This orthogonality condition 
cannot be detected point-wise in the energy density, but manifests itself as "miraculous 
cancellations" when the total energy is considered. In other words, to detect and account 
for this orthogonality requires a spatial mode-decomposition of the fields, as seen in [TQl 



That dominant energy condition cannot guarantee hyperbolicity can be seen also 
as a manifestation of this internal structure. Proposition [13] effectively says that 
perturbations ip which are parallel on N to the background solution <ft can also be 
controlled by the dominant energy condition. In the case where the dimension of N 
is greater than 1, it is precisely those additional degrees of freedom, which cannot be 
accounted for merely by the dominant energy condition, that requires the framework of 
regular hyperbolicity. 

5. Hyperbolicity of the Skyrme model 

It is easy to see that the Skyrme model can be considered as a theory of elasticity with the 
restrictions on the particle world-lines relaxed. This connection has also been observed 
by Slobodeanu |59j . From the point of view of relativistic fluids and elasticity, it then is 
perhaps less surprising that highly boosted Skyrmions are expected to be unstable [H] . 
Here we compute exactly the symbol m AB associated to the Skyrme model and exhibit 
an ultrahyperbolic type breakdown of hyperbolicity in tachyonic regimes. 

For the following discussion, we choose the scale factors to normalize the Skyrme 
Lagrangian [57JEE] to 



We ignore the mass term s(</>) as it plays no role in the discussion for hyperbolicity; 
here : M 1+3 — > § 3 is assumed as usual. Applying Theorem [6] we see that (with or 
without the positive mass term) the dominant energy condition is always satisfied for 
this model. 

Using (fl9l) . the coefficients m A b B can be computed to effectively be 



m% = g ab h AB (l + \d<j>\l h ) + d^Ad^B - g ah g cd d c <\> A d d <\> B - h A Bh CD d c <f d d <\>° (25) 



From Theorem [9] we already see that any time function of Minkowski space is an 
admissible time function. Therefore it suffices to examine the conditions for m^rja-rib 
to be positive definite. 

The following computations are most illustrative in a frame adapted to dcf). There 
are two cases: the kernel of dcf) being a degenerate (null) subspace, and everything else. In 
the everything else category, there exists an orthonormal frame of that simultaneously 
diagonalizes g and <p*h; in the case the kernel of dcj) is a degenerate, there exists an 
exceptional null frame (see, e.g. [27J Ch. 9, exercises 18 and 19] or [36], with some 



Ch. 5]. 
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cases ruled out by virtue of <p*h being positive semi definite) . The null case can be 
checked to be hyperbolic much in the same way as described below, so we shall omit 
its discussion (note that for the null case the canonical stress corresponding to <r 2 is 
degenerate; already it requires the presence of the wave-map term <j\ in the Lagrangian 
for hyperbolicity) and focus on the generic case of an orthonormal basis. Let the vector 
basis be e , ei, e 2 , e 3 , orthonormal with respect to the metric g, and diagonal with respect 
to {(jfh). Let their corresponding covector basis be / , fi, / 2 , f 3 where /^(e*) = 8y. And 
let Ao,Af, A|, A| be the values for (0*/t)(ej, ei) respectively. Note at least one of them 
must vanish from rank considerations. Let e\ be a unit vector in the tangent space of the 
target manifold such that the push-forward d<fi(ei) = A,e^ for such that d0( e «) = 0, 
choose e\ such that one of the e\ vanishes, and the rest forms an orthonormal basis for 
h. Take // to be their corresponding covectors. 
Then we can write 

9 ab = (-e ® e + ei ® e x + e 2 ® e 2 + e 3 ® e 3 ) ab (26a) 

h AB = (ft ® ft + ft ® /i + ft ® £ + ^ ® (266) 

<9 a A = (Ao/o ® e' + Aa/a ® + A 2 / 2 ® e' 2 + A3/3 <8> 4)a (26c) 

Let us consider, in this basis, the component m^ B . (The computations are symmetric 
in {1,2,3}, and the component m^ B is already treated by Theorem |9j) A direct 
computation shows that 

mf B = (1 + A 2 + A 2 )/o ® /o + (1 - Aq + Al)/i ® 

+(1 - A 2 + Ai)/ 2 ® ft + (1 - A 2 + A? + A 2 )^3 ® £ (27a) 

™ab = -A 3 A (/Da(/o)b (276) 

<b = A 3 Ai(/^)a(/0b (27c) 

Now first consider the case where Ao = 0. Then we see that m^ B is automatically 
positive definite. This case we have already discussed in the previous section: if d<p has 
a time-like element in its kernel, then that element is a valid observer field. 

For the case where Ao 7^ 0, we consider first when A3 = 0. Then we see that if 
Aq > 1 + min(A^, A 2 ), we have m^ B no longer positive. Furthermore, we have by the 
above computation, that in this case, for r\ = s/ 3 + rf , m a ^ B = s 2 m^ B + r 2 m^ B , a sum 
of a indefinite bilinear form with a negative definite one. Therefore we see that for any 
covector rj in the span of ft and ft, we cannot have m^BVaVb be positive definite. In 
particular, this implies that in this case regular hyperbolicity must fail. 

For the case where A ^ 7^ A 3 , and Ai = 0, we see that m^ B loses positivity as 
soon as A 2 , > 1. Now we consider again r\ = s/3 + rft. Computing explicitly we get 

m^VaVt = (s 2 - r 2 )(l + A 2 )/o ® ft + [s 2 (l -\ 2 + A 2 ) - r 2 (l + A 2 + A 2 )] ft ® ft 
+ [s 2 (l - A 2 ) - r 2 (l + A 2 )] ft ® ft + {s 2 - r 2 )(l + \l)ft ® ft 
-(s\ ft + r\ 3 ft)®(s\ ft + r\ 3 ft) 

where we see that the coefficient of ft ® ft will always be negative. This implies that for 
any 77 of this form, when Ai = and A 2 , > 1, we have that rriABVaVb cannot be positive 
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definite. 

From these computations, we get the following theorem 

Theorem 14. The Skyrme model is regularly hyperbolic when any of the following is 
true: 

(i) The kernel of dtp contains a time-like vector. 

(ii) The kernel of dip is a degenerate subspace. 

(Hi) The time-like eigenvector for has a corresponding eigenvalue with norm less 
than 1. 

The Skyrme model has an ultrahyperbolic breakdown of hyperbolicity when admits a 
time-like eigenvector with eigenvalue with norm greater than 1. 

Remark 15. The number "1" appearing in the third condition depends on the 
normalisation chosen for the Lagrangian, which dictates the interaction strength between 
the wave map and fluid-like terms. This result should be compared with the conclusion 
drawn by Crutchfield and Bell in [14], where their linear analysis shows that a sufficient 
criterion for break down of hyperbolicity is, in the above language, Aq > 1 + A^ + A2 + A3. 
So the present result sharpens the regime for which hyperbolicity fails. Note that this 
failure of hyperbolicity is not merely that of the regular hyperbolicity method. In parts 
of this regime it can be shown explicitly that the associated linear system is ill-posed; 
see the next remark. 

Remark 16. Fixing A 3 = and X\ < Aq — 1 < \\ (or Aq > 1 + \\ + A|, which is 
the Crutchfield-Bell condition), it is easy to check (using Descartes' rule of signs) that 
the linear, constant coefficient equation m^d^ip 3 = cannot be hyperbolic (in the 
sense that its polynomial symbol has the requisite number of real roots). Recall that 
hyperbolicity requires there to be a hyperbolic direction rj such that, for any £ transverse 
to r], the polynomial M(s) = det [m^^C + sr], ( + sr/)] has only real roots [T9U25] . We fix 
( = f 3 . Using that A3 = 0, we can exploit a symmetry condition that if 1] is hyperbolic, 
so will 7] — 2g(i], fs)f3. Using that the hyperbolic directions form a convex cone [19j, we 
can assume the would-be hyperbolic direction is orthogonal to ^3. For such an r], we 
have that M(s) is an even polynomial. The computations given before the statement 
of the theorem implies that lim s ^± 00 M(s)/ s 6 < 0, with M(0) < 0. Therefore M(s), a 
sixth degree polynomial, can have at most 4 real roots. The same argument, reversing 
the role of i], (, can also be used to rule out f$ as a hyperbolic direction, proving the 
claim. 

If we assume the local well-posedness result claimed in [10J, then the above theorem 
implies the following well-posedness property for the Cauchy problem of the Skyrme 
model. 

Corallary 17. The Cauchy problem for the Skyrme model with almost stationary initial 



data, where almost stationary is read to satisfy the hypotheses of Theorem^T^ is locally 
well-posed. In particular, a small perturbation of a static Skyrmion configuration gives 
rise to a well-defined evolution. 
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Remark 18. The existence of static configurations to the Skyrme model is a partially 
open problem. In certain symmetry classes the existence is known, see [TTII301 140J. In 
the case of the perturbation of a static configuration, local well-posedness also follows, 
by the computations above, using the techniques of [26]. The framework of regular 
hyperbolicity is not necessary in that regime. 

Remark 19. It has been pointed out to the author by Dan Geba that one automatically 
has local well-posedness with arbitrary smooth initial data, if one were to consider 
the spherically symmetric Skyrme model. Indeed, once the symmetry is imposed, the 
target space (being a quotient of a three dimensional manifold by a symmetry with 
two dimensional orbits) is effectively one dimensional, and the ultrahyperbolic-type 
breakdown which is due to the internal structure of the target manifold cannot occur 
(see also Section H~5]) . 

However, in this situation one has a hyperbolic illustration of the failure of the 
Coleman principle, analogous to that observed by Kapitanski and Ladyzhenskaya [30J. 
More precisely, the existence of uniqueness of solutions in the spherically symmetric 
class does not automatically guarantee that said solutions are in fact unique in the 
non-spherically symmetric class. This is because the spherically symmetric solution 
is only guaranteed to be stable under spherically symmetric perturbations. In our 
case, however, the above analysis shows that the linearized system is not necessarily 
stable under asymmetric perturbations, and therefore we cannot use Cauchy stability 
to conclude that symmetric initial data must lead to symmetric solutions! 
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